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If the energy density of the Universe before nucleosynthesis is dominated by a scalar field that 
decays and reheats the plasma to a temperature Trh smaller than the standard neutralino freeze 
out temperature, the neutrahno relic density differs from its standard value. In this case, the relic 
density depends on two additional parameters: Trh, and the number of neutralinos produced per (f) 
decay per unit mass of the (j) field. In this paper, we numerically study the neutralino relic density 
as a function of these reheating parameters within minimal supersymmetric standard models and 
show that the dark matter constraint can almost always be satisfied. 
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O ' I. MOTIVATION 

(N ■ 

The neutralino is considered a good dark matter candidate, one that naturally yields the required relic density. 
, ^ ■ Recently, however, it has been recognized that, owing to the experimental constraints and to the increased precision 
in the determination of the dark matter content of the Universe, the agreement between the observed relic density, 
^ ■ flcdmh^ = 0.113 ±0.009 and the relic density predicted with standard cosmological assumptions, is far from being 
[ a generic feature of supersymmetric models. In fact, models with bino-like neutralinos tend to overproduce them, 
I . and special mechanisms such as coannihilations or resonant annihilations are required to suppress the relic density 
^ ' down to the observed range. In contrast, models with higgsino- or wino-like neutralinos usually give too small a relic 
. density and, to compensate for it, large neutralino masses (m^ > 1 TeV) are needed. Non-standard cosmologies and 

■ in particular models with low reheating temperatures provide a plausible solution to this problem. These include 
', models with moduli decay 0, Q-ball decay 0j, and thermal inflation |^. In all of these models there is a late episode 

of entropy production and non-thermal production of the LSP in particle decays is possible. 

We concentrate on cosmological models in which the early Universe is dominated by the energy density of a scalar 
^) i field that after some time decays giving rise to the radiation dominated era. The decay of the scalar field into light 
' degrees of freedom and their subsequent thermalization -the reheating process- leaves the Universe at a temperature 
TjiH known as the reheating temperature. If, as assumed in the standard scenario, Tfuj is larger than the neutralino 
7*^' freeze out temperature (Tf.o. — the neutralino relic density is insensitive to its value. But, because we have 

Oh, no physical evidence of the radiation dominated Universe before big-baiig nucleosynthesis, Trh should be considered 
^ ' as a cosmological parameter that can take any value above a few MeV [a, Q ■ 

The existence of a weakly coupled scalar field that dominates the Universe during the process of neutralino pro- 
duction and freeze out may affect the relic density in several ways. It modifies the temperature-scale factor and 
the temperature-expansion rate relations j^, 0, ITol | that determine the freeze out condition. It dilutes the neutralino 
thermal abundance by increasing the entropy of the Universe with its decay into radiation. It may also increase the 

■ neutralino relic density by decaying into supersymmetric particles which in turn decay into neutralinos. As a result, 
the neutralino relic density may lie above or below its standard value depending on the reheating parameters. 

Several papers dealin g \y ith the neutralino relic density in the presence of a decaying scalar field already exist in 
the fitcrature US |lD,EllIllISIIlllllIEIlllIl- None of these, however, presents a systematic study of the 
combined effect of low reheating temperatures and non-thermal production, nor their implications within general 
supersymmetric models. The systematic study of thermal and non-thermal production mechanisms was introduced 
in Ref. Here we study the implications within the Minimal Supersymmetric Standard Model (MSSM). After 

reviewing the standard scenario, the equations that describe the reheating process will be introduced and numerically 
solved for different sets of parameters. We will elucidate the differences with respect to the standard computation and 
will investigate the effects of the reheating parameters on the neutralino relic density. Then, wc will compute Q.-^h'^ 
in generic MSSM supersymmetric models for different sets of the reheating parameters. Finally, a brief discussion of 
our results as well as our conclusions will be presented. 
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^ There is no qualitative change in our conclusions if we would use the recent more precise value f^cdrn'^^ = O-lO^io oo6' obtained for a 
ACDM model with scale-invariantprimordial perturbation spectrum through a global fit of cosmic microwave background, supernovae, 
and large scale structure data in IJ- 
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II. THE STANDARD SCENARIO 



In the standard scenario, it is assumed that neutrahnos are produced by scattering from the thermal plasma in a 
radiation dominated Universe. The neutralino relic density is then determined by the Boltzmann equation for the 
neutralino number density (n) and the law of entropy conservation: 

^ = -3Hn-{av)in^~nl^), (1) 



ds 
Uti 



-iHs . (2) 



Here t is time, s is the entropy density, H is the Hubble parameter, ri pn is the neutralino equilibrium number density, 
and {(Tv) is the thermally averaged total annihilation cross section j20j. Both the expansion of the Universe and 
the change in number density due to annihilations and inverse annihilations are taken into account in Eq. It is 
convenient to combine these two equations into a single one for Y = n/s and to use x = m/T, with T the photon 
temperature, as the independent variable instead of time: 

According to the Friedman equation, the Hubble parameter is determined by the mass-energy density p as 

3M 



H'^tttTTP^ (4) 



where Mp = 1.22 x 10"'^^ GeV is the Planck mass. The energy and entropy densities are related to the photon 
temperature by the equations 

P^^gMT)T\ s^^hMT)T' (5) 

where ^cff and hcs are the effective degrees of freedom for the energy density and entropy density respectively. If we 
define the degrees of freedom parameter gl^'^ as 

1/2 h^.ff ( . ^ T dh^ 



9oS 



then Eq. © can be written in the following way. 



dY / 45 \ gl^^m 



dx V 7rM|, 



az;)(r^-y4). (7) 



This single equation is then numerically solved with the initial condition Y = Yeq at a; ~ 1 to obtain the present 
neutralino abundance Yq. From it, the neutralino relic density can be computed 

atd/i' = ^ = "'^'°o^"^' = 2.755 X 10« Y^mJGeV , (8) 

Pc Pc 

where and sg Eire the present critical density and entropy density respectively. In obtaining the numerical value in 
Eq. (jSJ we used Tq = 2.726 K for the present background radiation temperature and /iofr(2o) = 3.91 corresponding to 
photons and three species of relativistic neutrinos. 

The numerical solution of Eq. jjjl shows that at high temperatures Y closely tracks its equilibrium value Y^q. In 
fact, the interaction rate of neutrahnos is strong enough to keep them in thermal and chemical equilibrium with the 
plasma. But as the temperature decreases, Ygg becomes exponentially suppressed and Y is no longer able to track its 
equilibrium value. At the freeze out temperature (Tf.o.), when the neutralino annihilation rate becomes of the order 
of the Hubble expansion rate, neutralino production becomes negligible and the neutralino abundance per comoving 
volume reaches its final value. 

From this discussion follows that the freeze out temperature plays a prominent role in determining the neutralino 
relic density. In the standard cosmological scenario, the neutralino freeze out temperature is about m^/20. In general, 
however, the freeze out temperature depends not only on the mass and interactions of the neutralino but also, through 
the Hubble parameter, on the content of the Universe. 
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III. LOW REHEATING TEMPERATURE SCENARIOS 



If the freeze out temperature is larger than the reheating temperature (Tf .o. > Tj^h) or equivalently if neutrahnos 
decoupled from the plasma before the end of the reheating process, the neutralino relic density will differ from its 
standard value. In that case, the dynamics of the Universe is described by the following equations 



— = -3Hp^-T^p^ (9) 
^ = -3Hn - (av) {v? - n%) + ^^^P4> (10) 



|.-3i.. + ^ (11) 

where m^, T^, and p^ arc respectively the mass, the decay width and the energy density of the scalar field, and h 
is the average number of neutralinos produced per (f> decay. Notice that b and enter into these equations only 
through the ratio h/m^ and not separately. Finally, the Hubble parameter, H , receives contributions from the scalar 
field. Standard Model particles, and supersymmetric particles, 

^^^{p^ + PSM+Px)- (12) 

These equations take into account the evolution of the scalar field, its contribution to the entropy density, and its 
possible decay into supersymmetric particles. As expected, in the limit p^ ^ Q they reduce to the standard scenario. 
In writing them, we have implicitly assumed that initially the 4> field inflated the Universe to a state of negligible 
entropy, that in each (p decay an energy is contributed to the matter and radiation plasma, and that the decay 
products rapidly thermalize (i.e. reach kinetic, although not necessarily chemical, equilibrium). 

In cosmologies with a decaying scalar, therefore, the neutralino relic density depends on two additional parameters: 
h/ms, and F^. In specific models describing the physics of the (/> field, these parameters may be calculable (see 
Ref. [13 and references therein). Our approach, however, will be completely phenomenological. Wc will treat fo/m^ 
and F^ as free parameters subject to the constraints that the dark matter bound imposes on them. For convenience 
instead of b/m^ we will sometimes use the dimcnsionless quantity 

, /lOO TeV\ 

V^bi (13) 

\ J 

as the free parameter. Clearly, b = t] for = 100 TcV. F^, on the other hand, is commonly expressed in terms of the 
parameter Thh ~ conventionally called the reheating temperature - defined by assuming an instantaneous conversion 
of the scalar field energy density into radiation. 



The Eqs.|nj^2|above assume thermal equilibrium. As explained in Ref. [l9l |. this assumption is well justified. Let us 
see why. First note that during the epoch in which the Universe is dominated by the decaying field, H is proportional 
to 0. We can derive this dependence starting from the evolution equation of the entropy per comoving volume 
S = sa^, which is eq^Jniultiplied by a'^, namely dS/dt ~ T^p^a^ /T . During the oscillating </> dominated period, pcj,a^ 
is constant, thus a? oc p~^^ . Using H oc ^/p^ (xt~^ we obtain oc and writing T (x t" we get S ~ T^a^ oc 
Substituting these expressions for S and T as functions of time in the evolution equation for S, and matching the 
powers of t on both sides, determines a = —(1/4). Hence, during the oscillating cj) dominated epoch H oc t^^ oc T"*, 
and p0 oc iJ2 oc T^. Since H ~ T^^/Mp at T = Trh, it follows that H ~ T'^ / {T'^^Mp). 

Neutralinos are in thermal (also called kinetic) equilibrium with the radiation when their scattering rate off rela- 
tivistic particles is faster than the Hubble expansion rate, F'^'^^" > H. Because F'^'^'^** ~ 'T^70'scatt — T^cfscatti we see 
that kinetic equilibrium is maintained at low temperatures, T < r^jjA/pCTgcatt- This is contrary to the usual radiation- 
dominated scenario in which the temperature has to be large enough for kinetic equilibrium to be maintained. The 
origin of this difference is the strong dependence of H on the radiation temperature. Using the known relation 



av) 



Q.stdh — 10 GeV / {av), and taking Cscatt of the same order of magnitude as the annihilation cross section 
gives T < (10^ MeV /nstdh^){TuH/McV)^ for kinetic equilibrium to be maintained. For the vaUdity of Eas.lMTH we 
need to assume kinetic equilibrium before neutralino production ceases: the right hand side is larger than Trh for 
^stdh^ ^ 10^(TRH/MeV), so wc safely assume that kinetic equilibrium is reached. 
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To solve the system of equations (|^- H12() . besides rj and Tuh we have to specify the initial conditions for p^, n^, 
and s. Giving the value Hj of the Hubble parameter at the beginning of the (j) dominated epoch amounts to giving the 
initial energy density p^j in the (j) field, or equivalently the maximum temperature of the radiation Tmax- Indeed, one 
has Hj ~ plj/j/Mp ~ T^^-^/ (T^^Mp). The latter relation can be derived from ~ ^^/^rh ^^'^ ^^"^ consideration 
that the maximum energy in the radiation equals the initial (maximum) energy p^,/. At the fundamental level the 
initial energy density of the scalar field depends on the particular model describing it. In hybrid models of inflation, 
for instance, p^j ~ m^, whereas a much larger value is expected in chaotic inflation p^j ~ m^Mp. From a physical 
point of view it is not surprising that the relic density is sensitive to the initial (f) energy density, for it determines 
the available energy of the Universe. If p^j is too small, the Universe can never get hot enough to thermally produce 
neutralinos and the relic density is typically suppressed. If the neutralino reaches chemical equilibrium, it is clear that 
its final density does not depend of the initial conditions. An approximate condition for reaching chemical equilibrium 
is ^3 i'^^) ^ 10~^GeV~^ (TOjj,/100GeV)(rRH/McV)^'^. Even without reaching chemical equilibrium, the neutralino 
density is insensitive to the initial conditions provided the maximum temperature of the radiation Tmax ^ "m-x C3 ■ 

Regarding the other two initial conditions, we assume that at an arbitrary reference temperature Ti the initial 
entropy is given by Eq. and that the neutralino number density is negligible, rii = 0. We explicitly checked that 
when chemical equilibrium is reached, using different values for T, and rii does not modify the neutralino relic density. 
Thus, our results will all be independent on the initial conditions. 

At early times, when tV^ <^ 1 and the 4> field is oscillating around the minimum of its potential, the scalar field 
energy density per comoving volume p^a^ is essentially constant and useful analytical results can be derived from eqns 
(P))- (|12|l . The temperature and the scale factor are related by T oc a^"^/® whereas the temperature and the expansion 
rate are linked hy H oc T^, marking the departure from the standard cosmological scenario -where T oc and 
H ocT^. This is not the only relevant departure, though. From a qualitative point of view, at least three important 
differences between the standard cosmology and the scenario with a decaying scalar field can be identified. They are 

• During freeze out the Universe may not be dominated by radiation. Initially, the energy density is dominated 
by the scalar field and only at the end of the reheating process the Universe becomes radiation dominated. 

• The entropy density of the Universe is not conserved but increases during the reheating era due to the scalar 
field decay. 

• Neutralinos could be produced non-thermally in cf) decays. This production mechanism could even be the 
dominant source of neutralinos. 

In Ref. ^3 the solutions to Eqs. (Pl <12|l were obtained analytically (as well as numerically). There are four different 
cases which result from the different ways in which the density fl^h^ depends on Trh: 

Case (1), thermal production without chemical equilibrium. In this case il^ oc T^jj- relic density was estimated 
in Ref. 



^ / TT,,^ \ ^ f 0^ ^^"^ 



»x / 100GeV \ VTrhV _ 



(15) 



Case (2), thermal production with chemical equilibrium. In this case il^ oc T^jj. The neutralino freezes out while 
the universe is dominated by the (f> field. Its freeze-out density is larger than usual, but it is diluted by entropy 
production from <f> decays. The new freeze-out temperature T^^^ is determined by solving njav) ~ if at T = T^^^ . 
Using the relations between H, a, and T in the decaying-(/) dominated Universe, one finds [3. Il0|| 

o^c.r3^rf.o.(Tf^EW)-4^^^^. (16) 

Our numerical results indicate a slope closer to T^jj. at least partially due to the change in T/^^w (g^^^ pjg^ |21below) 
Case (3), non-thermal production without chemical equilibrium. Here oc Trh. Non-thermal production is not 
compensated by annihilation. The production of neutralinos is purely non-thermal and the relic density depends on rj. 
It can be estimated analytically as follows. For each superpartner produced, at least one LSP will remain at the end 
of a chain of decays (due to i?-parity conservation), and thus ~ bn^. Here = p^/m^. At the time of 0-decay 
~ m^bp^/m^ ~ ^RH' ''^^'^ t^^6 entropy is s ~ T^^h- Hence p^/s ~ Trh and Yq = Idecay — bT^n/rii^. It follows 
that, 

o ir,3„ f "^x A f ^RH 



2 X lO'^Ty i — (17) 

..cdm VlOO GeV/ \McVj ^ ' 

Case (4), non-thermal production with chemical equilibrium. In this case f2 oc Tj^jJ. Annihilation compensates for 
the non-thermal production of neutralinos until the non-thermal production ceases at T = Trh. The condition for 
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determining the relic density is T™" ~ at T = Trr. This leads to Fq ~ Irh ^ T^/{sRH{crv)) ~ l/(TRHMp(cru)). 
From here it follows that 

~ (Tf.o./rRH)f^std. (18) 

In R,ef . [T^ . it was concluded that only neutralinos whose standard relic density is fistd ^ 10~^(100GeV/TO^) cannot 
be brought to f2cdmj independently of 77, while all neutralinos with larger f2std can be made to account for all of cold 
dark matter within these scenarios. In order to have an idea of how many neutralino models cannot be brought to 
have the full dark matter density, Fig. ^ shows the range of values of figtd for a few million mSUGRA models (data 
courtesy of Ted Baltz). One can see that very few of them, namely only some with ~ mz/2 for which neutralinos 
annihilate resonantly through a Z boson, have ilstd ^ 10~^(100GeV/m^), i.e. are below the slanted line in the figures. 
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FIG. 1: Range of values of the standard relic density for mSUGRA models. Only those with Qstd ^ 10-^(100GeV/m^), i.e. 
those below the slanted line, cannot be brought to have the total density of dark matter in low reheating temperature scenarios. 

In the following we investigate quantitatively how the reheating process modifies the evolution of the neutralino 
abundance during the reheating era. To numerically integrate the equations describing the reheating process, we 
modified the DarkSUSY code so that the three equations are solved simultaneously. The advantage of 

using a program like DarkSUSY is twofold. On the one hand, it provides an efficient and precise algorithm to solve 
differential equations like Hl()|l . On the other hand, it automatically computes the total annihilation cross section (av) 
in terms of supersymmetric parameters. 

To begin our analysis, let us consider thermal production only {rj = 0) and let us see how low reheating temperatures 
affect the neutralino relic density. In Fig. |21we plot the evolution of the neutralino density Y = n/s as a function 
of T for a particular choice of mSUGRA parameters: M1/2 = rriQ = 600 GeV, Aq = 0, tan/3 = 10, and /i > 0. 
For this choice, flstdh^ = 3.6, thus we are examining the case of a slightly overdense neutralino. The solid lines 
show the equilibrium density and the neutralino density in the standard scenario. The other lines show the effect of 
different reheating temperatures between 10 GeV and 10 MeV. Note that the value of the parameter Trh, shown by 
a short vertical line in each curve, coincides well with the actual reheating temperature (at which the curve changes 
slope, indicating that the oscillating (p dominated era finishes and the radiation dominated era begins). We are in 
our case (2), thermal production with chemical equilibrium, for which f2^/i^ cx T^^j, as verified in the figure. The 
curve for Tjih = 10 GeV closely tracks the standard one, yielding an abundance just below the standard prediction. 
As we move to smaller reheating temperatures, the freeze out occurs earlier, at a temperature r^EW ^ Tf a. and so 
the neutralino abundance at freeze out is larger than in the standard scenario. But, between T^'^"^ and Trh the 
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FIG. 2: The evolution of the neutrahno abundance for different values of Trh and 77 = (i.e. with only thermal production). 
In this figure we use the mSUGRA parameters M1/2 = 'mo = 600 GeV, Aq = 0, tan/3 = 10, and fi> 0. The neutralino mass is 
= 246 GeV. The 




FIG. 3: The evolution of the neutralino abundance for Trh = 1 GeV and several values of 77. Minimal SUGRA parameters as 
in Fig. H 



neutralino abundance is diluted away by the entropy released during the decay of the scalar field. In this regime, 
represented by the descending straight lines, Y (x T^. Indeed, n^a^ (with a the scale factor) remains constant 
whereas from a oc follows that sa^ decreases as T~^; hence Y = n^/s oc T^. Because the entropy dilution 

effect is dominant, the neutralino relic abundance ends below the standard prediction and the smaller the reheating 
temperature, the smaller the relic density. Thus, if neutralinos are produced thermally, a low reheating temperature 
entails a suppression of the relic density with respect to its standard value. 

Non-thermal production of neutralinos, however, can compensate for such suppression. In Fig. |21 we show the 
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FIG. 4: The neutralino relic density as a function of 77 for several values of Tbh- Minimal SUGRA parameters as in fig|21 
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FIG. 5; The neutralino relic density as a function of Tb,h for different values of 77. Minimal SUGRA parameters as in fig|5| 



evolution of the neutralino abundance for Tj^h = 1 GeV and different values of 77. The dependence of Y on the 
temperature between Tf .o. and Te.h is now different due to the production of neutralinos in the </> decay. As we go 
from to larger values of 77, the neutralino abundance increases and at 77 ~ 10~^ it becomes larger than the standard 
prediction. We are in our case (3), non-thermal production without chemical equilibrium, in which production is not 
compensated by annihilation. In Eq. E|we see that for fixed Thh the final relic abundance increases with 77, Vl^ cx 77. 
Setting — r^cdm, we find the critical value of 77 above which the relic abundance is larger than its standard value. 
Replacing Trh = 1 GeV, mx = 246 GeV and ficdm/i^ = 0.11, we find that ^-^Ji^ becomes larger than ^stdh^ — 3.6 for 
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FIG. 6: Regions in the plane {TRH,r]) compatible with the WMAP range for two different supersymmetric spectra. The solid 
region corresponds to mo = M-1/2 = 400 GeV, tan/3 = 10, Aq = 0, and /i > 0, parameters which yield Q.^h? = 1.8 and = 160 
GeV . The hatched region corresponds to mo = 600 GeV, Mi/ 2 = 800 GeV, Ao = 0, tan/3 = 45, and ^ < 0, parameters which 
yield Q.^h^ = 0.02 and m^ = 333 GeV. 



the critical value rjc ~ 0.7 x lO""'. For larger values of 77, the density of neutralinos is large enough for annihilations 
to become important. We are then in our ease (4), non-thermal production with chemical equilibrium, and Eg. 1181 
gives the final abundance. Equating Eos . II 71 and 1181 gives the value of -q at the cross-over, which is 



This cross-over value of 77 is about 0.8 xlO~^ for the example we are considering. Thus the curve with = 10~^ 
in Fig. 13 corresponds to case (4). For large values of ry, the non-thermal production of neutralinos is efficient and 
Y cx i-jT for T > Tbh (as follows from = n^a^ cx -qt cx i]T^'^). For rj = 10~^, for instance, the relic density is more 
than two orders of magnitude larger than the standard one. Hence, if neutralinos are produced non-thermally a low 
reheating temperature can yield a relic density below or above the standard result, depending on the values of Trh 
and rj. 

What can be constrained with the observations today, however, is not the evolution of the neutralino abundance 
but just its asymptotic value -the neutralino relic density. In Fig. 0] we show fl^h^ as a function of ?y for different 
values of Trh- In the figure, Tf.o. — 12 GeV. For T^h = 50 GeV (solid fine), the rehc density is independent of b and 
equal to its standard value indicating that the freeze out took place after the reheating era, in a radiation dominated 
Universe. For smaller values of Tf^j, however, fi^/i^ does depend on rj. This is the regime of case (3), Eq.EI in which 
non-thermal production is not compensated by annihilation. As mentioned above, setting = l^cdm in Eg. 1171 one 
finds for each value of T^h the critical value of 77 above which the relic abundance is larger than its standard value. 
For the particular model in the figure this value is rjc ~ 0.7 x 10~^(MeV/TR^). Eg. 1191 gives the cross-over value 
T]c-o — 0.8 X 10^ {MeV /Trh)^ at which annihilation becomes important, and for larger values of rj the relic density is 
given by Eg.lTHlfcase (4)). Notice that rjc ~ ??c-o for Trh — lOGeV. 

It is clear from Fig. ^ that the relic density is below the standard result for small 77, it increases with ry and 
at a certain point, 77 = rjc, coincides with the standard one; for rj > 77c, the relic density is larger than in the usual 
cosmology. The value of 77c is inversely proportional to T^h at low reheating temperatures, Tfin < IGeV in the figure. 
In this regime, the relic density becomes a straight line signaling dominant non-thermal production of neutralinos 
without chemical equilibrium (case 3, Eg. 1171 f2^ft,^ oc 77). 

In Fig. |S1 the relic density is shown as a function of T^/f for different values of rj. This figure is similar to Fig. 1 of 
Ref. [13 ■ We observe that at small Trh, Tjih < IGeV in the figure, the relic density is proportional to Tj^h- This is 
the regime of non-thermal production without chemical equilibrium, in which fl^h^ cx rjTun (case (3), Eg. I17II . As 
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FIG. 7: Scatter plot of the neutralino relic density as a function of the neutralino mass in the standard scenario. 

TjiH increases above IGeV there are in the figure some descending and some ascending curves. The descending curves 
correspond to non-thermal production with chemical equilibrium, in which fl^h^ oc I/Tr// (case (4), Eg. I18() . The 
ascending curves correspond to thermal production with chemical equilibrium, where n^h'^ oc T^h (case (2), Ea. [T?)|l . 
At TjiH > Tf.o. — 12 GcV in the figure, fl^h^ is independent of and Thh, and equal to its standard value. As a 
general conclusion, depending on the particular choices for Thh and rj, the relic density may be larger or smaller than 
its standard value. 

Though specific values for the supersymmetric parameters were chosen in Figs.0El tlicy actually represent a generic 
situation whenever l^std is small enough for chemical equilibrium to be reached at high Tf/n- Taking different points 
in the parameter space is equivalent to rigidly translating the figures without altering their shapes. 

Another way of visualizing the effects of Tf/n and rj on the relic density is by studying the regions compatible with 
the WMAP range, ficdm^^ = 0.09 — 0.13. In Fig.Elwc show such regions in the plane {Tj^h, v) for two different sets of 
mSUGRA parameters. Neither model is viable in the standard cosmological scenario because one overproduces dark 
matter (solid region, corresponding to flstdh'^ = 1-8, = 160 GcV) and the other underproduces it (hatched region, 
corresponding to fistd^^ = 0.02, — 333 GeV). But, as this figure illustrates, in cosmologies with low reheating 
temperatures the neutralinos in both models can easily account for the dark matter content of the Universe. Notice that 
the first model requires a reheating temperature below 4 GeV. This is the value of Trh above which thermal production 
dominates and the relic density becomes independent of rj and proportional to T^jj (corresponding to our case (2)) 
until Trh reaches the standard freeze-out temperature Tf.o. — rfiy^jlQ ~ 8 GeV, at which point the density becomes 
the standard rehc density. Thus here (0.11/ristd^^)^^'^Tf.o. — Tf.o./2 ~ 4 GeV. The second model demands a value of 77 
above lO"'^. It is also evident from the figure that in underabundant models for each value of 77 two different reheating 
temperatures can be found that give the correct relic density. In fact, in Ref. [l9l | it was shown that for neutralinos with 
standard densities ficdm ^ f^std ^ 10~^(100GeV/TO^), there is no solution for 77 < 10-^(100GeV/TOx)2(17^dm/f^std), 
there are two solutions for 10~^(100GcV/7Ti^)^(ricdm/^^std) ^ V ^ 10^^(100GcV/to^); and there is a single solution 
for larger values of 77. 

IV. THE NEUTRALINO RELIC DENSITY IN SUPERSYMMETRIC MODELS 

In the MSSM, neutralinos are linear combinations of the fermionic partners of the neutral electroweak bosons, called 
bino (i3°) and wino (1^3), and of the fermionic partners of the neutral Higgs bosons, called higgsinos (-ff°, H^)- We 
assume that the lightest neutralino, x, is the dark matter candidate. Its composition can be parametrized as 
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FIG. 8: Scatter plots of the relic density as a function of the neutralino mass for i] — and different values of Trh- 
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FIG. 9: As fig 8 but 77 = 10"^ 

Because the neutralino interactions are determined by its gauge content, it is useful to distinguish between bino-like, 
wino-likc, and higgsino-likc ncutralinos according to the dominant term in H20|l . 

Bino-likc ncutralinos annihilate mainly into fcrmion-antifcrmion pairs through sfcrmion exchange, thus they typi- 
cally have a small annihilation cross-section that translates into a standard relic density larger than observed. Agree- 
ment with the observed dark matter abundance can still be achieved in standard cosmological scenarios but only in 
restricted regions of the parameter space where special mechanisms such as coannihilations or resonant annihilations 
help reduce the relic density. Bino-like ncutralinos are a generic prediction of minimal supergravity models. 

Wino-like and higgsino-like ncutralinos, on the other hand, annihilate mostly into gauge bosons (VF+VK^, ZZ, if 
kinematically allowed) through neutralino or chargino exchange; otherwise they annihilate into fermions. Through 
coannihilations with ncutralinos and charginos of similar mass, their standard relic density is rather small. Neutralino 
masses as large as 1 TeV for higgsinos or 2 TeV for winos are required to bring their thermal density within the 
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FIG. 10: Scatter plots of the relic density as a function of the neutralino mass for Trh = 1 GeV and different values of 77. 
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FIG. 11: As Fig. 10 but for Trh = 10 MeV. 

observed range as can be seen in Fig. [7| (see below) . Wino-likc and higgsino-like neutralinos can be obtained in models 
with non-universal gaugino masses; AMSB models, for instance, feature a wino-likc neutralino. 

Since it is our aim to study the effect of a decaying scalar on the neutralino relic density within supersymmetric 
models, we must consider models generic enough to allow for bino-like, wino-like, and higgsino-like neutralinos. To 
that end, we will examine MSSM models defined in terms of the parameter set M3, Af2, Mi, niA, ^J■, tan/?, toq. At, 
and Ai,. Here Mi are the three gaugino masses, is the mass of the pseudoscalar higgs boson, and tan/3 denotes the 
ratio V2/V1. The soft breaking scalar masses are defined through the simplifying ansatz Mq = Mjj = Md ^ Me ^ 
Ml = mo whereas the trilinear couplings are given by Ajj = diag(0, 0, At), Ad = diag(0, 0, Ab), and Ae = 0. All 
these parameters are defined at the weak scale. Specific realizations of supersymmetry breaking such as mSUGRA or 
AMSB are particular examples of these models. 
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We perform a random scan of such parameter space within the following ranges 

10GeV< Af„TOA,/i <50TeV (21) 

10 GeV < mo < 200 TeV (22) 

-3mo < A,Ab < 3too (23) 

1< tan/3 < 60 (24) 

A logarithmic distribution was used for Mi, niA, and mo, and a linear one for At, Ah, and tan/3; the sign of /i 
was randomly chosen. Accelerator constraints (as contained in DarkSUSY version 4.1 |2j|) were imposed on these 
models. For each viable model generated in this way, we compute the neutralino relic density for different values of 
the reheating parameters rj and Tjih (we computed 1,700 models for each pair of values of these two parameters and 
each model is a point in the scatter-plots). In this section we present and analyze such results. 

We start by showing in Fig. [7| the relic density as a function of the neutralino mass in the standard cosmological 
scenario. The horizontal band indicates the WMAP range. From this figure we can easily tell the different kinds 
of neutralinos: the relic density of wino-like and higgsino-like neutralinos is determined essentially by the neutralino 
mass because their annihilation cross section is insensitive to the sfermion and higgs sectors. Consequently, models 
featuring these neutralinos are distributed along straight lines and are clearly identifiable in the figure. The line 
corresponding to wino-like neutralinos crosses the WMAP interval at m^ ~ 2 TeV, whereas that corresponding to 
higgsino-like neutralinos does it at m^ ~ 1 TeV, because winos annihilate more efficiently than higgsinos. Since the 
bino annihilation cross section depends on the sfermion spectrum, the relic density in models with bino-likc neutralinos 
can vary over a large range even for a given neutralino mass. As seen in the figure, bino-like neutralinos typically 
yield too large a relic density in the standard cosmological scenario. 

Let us now consider models with a decaying scalar. A scatter plot of the relic density as a function of the neutralino 
mass for 77 = and different values of Trh is shown in Fig. |S| The suppression of the relic density due to small 
reheating temperatures is evident in this figure. For Trh = 10 GeV, this effect is noticeable for large neutralino 
masses; the wino and higgsino lines bend downward close to m^ ~ 300 GeV, and heavy binos are brought closer to 
the WMAP range. For Tjih = 1 GeV the suppression is larger and only light binos (m^ < 200 GeV) give a relic 
density compatible with the observations. If Tj^h < 100 MeV the neutralino relic density is too small to account for 
the dark matter of the Universe. 

In general, however, low reheating temperatures can be accompanied by non-thermal production of neutralinos, 
modifying the previous results. In Fig. |^ we show the relic density for 77 = 10~^ and different reheating temperatures. 
The straight bands with a linear dependence on m^ shown in this and subsequent figures corresponds to our case 
(3) solutions, non-thermal production without chemical equilibrium, given in Eq. 1171 For Tjih = 10 MeV neutralino 
production is entirely non-thermal without chemical equilibrium and consequently the relic density is proportional 
to the neutralino mass fEa. I17|l . Obtaining the correct relic density requires a neutralino mass of about 4 TeV. If 
Trh = 100 MeV thermal effects are relevant only for hght winos (m^ < 200 GeV, for which Tf.o. ~ 10 GeV < Trh); 
heavy winos as well as binos and higgsinos are produced non-thermally. The neutralino relic density lies within the 
WMAP range only for m^ ~ 400 GeV. For Trh = 1 GeV, light (m^ < 1 TeV) binos, winos and higgsinos become 
clearly distinguishable. The correct relic abundance can be obtained at m^ < 100 GeV for binos, m-^ ~ 250 GeV 
for higgsinos, and m^ ~ 450 GeV for winos. If m^ > 1 TeV neutralinos are produced non-thermally and have a 
relic abundance well above the WMAP range. Finally, if TfiH = 10 GeV attaining the observed dark matter density 
requires m^ < 300 GeV for binos, ~ 500 GeV for higgsinos, and m-^ ~ 800 GeV for winos. 

We will now look at the effect of different values of r] for a given Tub- In Fig- El we show the relic density as a 
function of the neutralino mass for Trh = 100 MeV and several values of 77. For 77 = 0, the relic density is suppressed 
lying well below the WMAP range. For 77 = 10^^, the WMAP density can be achieved outside the mass range 
explored numerically along the extrapolation of the slanted straight line to higher masses; the crossing occurs at 



m 



X 



5 X 10^ GeV, as follows using Eq. El With 77 = 10~^, the correct relic abundance can be achieved only for 
neutralino masses between 300 and 400 GeV. With 7/ > 10""^, only a wino-like neutralino with m^ ~ 100 — 200 GeV 
can account for the dark matter. Fig.^Jis as Fig. 10 but with Tjih = 10 MeV. For such a small reheating temperature, 
the thermal production of neutralinos is negligible. The dispersion of points discernible in this figure for 77 = 1/2 (and 
for 77 = 10^^) are due to the annihilation of neutralinos produced in the decay of the field. These points correspond 
to non-thermal production with chemical equilibrium, our case (4) , eq. 1181 It is still possible to achieve the observed 
relic density for 77 between 10^® and 10^^. 

To summarize graphically the effects of Thh and 77 on the relic density in generic supersymmetric models, we 
present in Fig. 1121 an array of scatter plots of the relic density for different reheating parameters. Each panel includes 
the same supersymmetric models and so the differences in the relic densities are entirely due to the effects of rj and 
Trh- The reheating temperature is constant along columns while 77 is constant along rows. Several conclusions can be 
drawn from this figure. The fourth row (7/ = 10""^) tells us that, if 7/ = 10"'^ bino-like neutralinos arc disfavored due 
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FIG. 13: r; as a function of the neutralino mass for models compatible with the WMAP range. 



to their large relic density, and the dark matter is most easily explained with wino or higgsino-like neutralinos. From 
the last column {Tuh = 10 MeV), we observe that if Tuh ~ 10 MeV, r/ should be between 10~^ and 10~^ in order 
to give the correct relic density. The crucial observation, however, is that by varying rj and T^h "we can move all the 
points above {Tj^h = 10, ?? = 1/2) or below (Trh = 100 MeV, 77 = 10^^) the WMAP range. For any supersymmetric 
model, therefore, it is possible to find values of rj and TfjH that give the observed value of the relic density. In other 
words, by choosing rj and T^jf appropriately the dark matter bound can always be satisfied. 

A question we have not directly addressed so far is what the parameter space compatible with the WMAP range 
is. Given the complexity of the parameter space involved (10 supersymmetric parameters plus Thh and b), it is 
not possible to answer this question in all its generality. What we can do is fix some parameters and see if there are 
constraints on any others. In Fig.^l we show the value of rj required to obtain 0.09 < ft^h^ < 0.13 as a function of the 
neutralino mass, for three different values of the reheating temperature, Tj^h = 1 GeV, 100 MeV, 10 MeV. Again here 
the straight bands proportional to 'm~^ are due to non-thermal production of neutralinos without chemical equilibrium 
(our case (3), eg. I17|) . For TfjH = 10 MeV, most neutralinos are produced non-thermally and the relic density is 
proportional to 77 fEa. I17|l . Heavy neutralinos {rriy^ > 1 TeV) require 77 < 4 x 10~^ whereas the lightest neutralinos 
demand i] ~ lO"**. From this band including most models, a small branch of wino-like neutralinos deviates (at 
rriy, ~ 100 GeV, 77 ~ 4 x 10~^) and reaches values of 77 as large as 10~^. These models have a very small standard relic 
density and reach the required dark matter density range through the solutions of our case (4), Eg. 1181 non-thermal 
production with chemical equilibrium, which require large values of 77. For Tfuj = 100 MeV, the figure is similar with 
the wino deviation occurring at larger masses (771^ ~ 250 GeV) and 77 extending up to 10"'^. For Tjih = 1 GeV, on 
the other hand, two different branches corresponding to winos and higgsinos can be distinguished. Moreover, some 
points lying below the non-thermal production band are also present. 



V. CONCLUSIONS 



In cosmologies with a decaying scalar field, the dark matter bound on supersymmetric models can almost always 
be satisfied. In fact, by adjusting Trj/ and rj the relic density can always be brought into the observed range except 
for neutralinos with ^Igtd ^ 10~^(100GeV/m^). 

Models that in the standard cosmological scenario have a relic density above the WMAP range, as is generically 
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FIG. 14: The lines show the values of mo and M1/2 for which the relic neutralino density has the central value of the range 
imposed by WMAP bounds, namely Qcdmh^ = 0.11 for different values of 77. Minimal SUGRA model with Ao = 0, tan/3 — 10 
and /i > 0. 

the case for bino-like neutralinos, can be rescued by decreasing the reheating temperature of the Universe (either just 
suppressing the thermal production, if Tnu is below but close to the standard freeze-out temperature, or producing 
them non-thcrmally with a non- unique combination of Tfin and rj for lower Tj^h values). Models that in the standard 
scenario have a relic density below the WMAP range, as for higgsino- and wino-likc neutralinos, can be saved with 
an appropriate and non-unique combination of Trh and rj. 

For dark matter detection purposes, the ideal neutralino - i.e. one that accounts for all of the dark matter content of 
the Universe, has a relatively small mass, so that its number density is large, and has large interaction and annihilation 
cross sections - is an oddity within the standard cosmological scenario. A light neutralino with a large interaction cross 
section typically yields a relic density smaller than the dark matter density. Thus, within the standard cosmological 
scenario, we are esentially led to two possibilities: a light neutralino with a small annihilation cross section, such as a 
bino-like neutralino; or a heavy neutralino (m^ > 1 TeV) with a large annihilation cross section, such as a higgsino or 
wino-like neutralino. If neutralinos are produced non-thermally, as is possible in low reheating temperature scenarios, 
the neutralino relic density is determined by the physics at the high scale, which determines Tjih and 77, and not by 
the annihilation cross section. Therefore, by adjusting Tjih and rj, it becomes possible for a light neutralino with a 
large annhihilation cross section to account entirely for the dark matter of the Universe. Ultra-heavy neutralinos, in 
the TeV mass range, become good dark matter candidates too. 

For example, the usual narrow corridors in toq, M1/2 space of good dark matter neutralino candidates in mSUGRA 
models are replaced by other narrow corridors, which, however, depend on the physics at the high scale that determines 
T]^H and r; = b/m^. The toq: -^1/2 plane for a mSUGRA model with Aq ~ 0, tan/3 = 10 and fi > is shown in 
in Fig. 1141 The gray (green) area is forbidden by experimental bound. The almost vertical lines in the figure show 
where neutralinos have ilcdm^^ = 0.11, the central value of the range imposed by WMAP bounds, for Trh — IGcV 
and different values of t]. We see that in all points in the toq, M1/2 space in Fig. 1141 neutralinos can account for all of 
the cold dark matter, with appropriate values of Trh and 77. In other words, in low reheating temperature scenarios 
with a decaying scalar field, cosmological data constrain models at the inflationary or Planck scale which determine 
Trh and 77. 

Future accelerator or dark matter detection experiments might find a neutralino in a region of the supersymmctric 
parameter space where its standard relic density is larger than the observed cold dark matter density. This would 
tell us that the history of the Universe before Big Bang Nucleosynthesis does not follow the standard cosmological 
assumptions, and would point towards non-standard scenarios, such as those with low reheating temperature studied 
here. 
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